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In this talk I discuss the antenna subtraction method for isolating infrared (IR) singularities of 
jet cross sections in perturbative QCD. The method is applied at next-to-next-to-leading order 
(NNLO) to dijet production in hadron collisions at the LHC. The double real radiative corrections 
to the dijet cross section are considered and their IR behaviour is examined. IR subtraction terms 
are constructed to absorb numerical divergences of the cross section in the single and double unre- 
solved regions of phase space using the antenna subtraction method. A pictorial representation of 
matrix elements and antenna functions is presented with specific examples of how such diagrams 
can be used in practical calculations. 
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1. Introduction 



One of the most powerful tools for describing the physics of high energy hadronic collisions is 
perturbative QCD in the context of the parton model. In this framework, the cross section for two 
partons to scatter into a number of jets can be calculated as a perturbative expansion in the strong 
coupling, a s [1], 
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where £1,^2 an d fi(^i),fj{^2) denote the momentum fraction and Parton Density Function (PDF) 
associated with initial state parton i and j respectively. By calculating the higher order contribu- 
tions, the theoretical prediction for the cross section can be improved in several respects [2]. In 
high energy collisions where a s . is small, higher order corrections naturally improve the conver- 
gence of the perturbative series and significantly reduce the renormalization and factorization scale 
uncertainty. The partonic cross section can only be compared with experimentally observed jets by 
use of a jet clustering algorithm. Such algorithms benefit from having more particles in the final 
state and so higher order calculations allow a more realistic comparison between parton level and 
hadron level jets [3]. Similarly higher order corrections allow for multiple initial state radiation 
which in turn generates transverse momentum for the final state. 

Next-to-leading order (NLO) calculations of the dijet cross section have been available for 
some time and have been tested against experimental data at both the TEVATRON and the LHC [4, 
5]. Currently the cutting edge for jet cross sections is the NNLO conuibution to dijet production 
with hadronic initial states. NNLO calculations contain several components: the two loop double- 
virtual cross section, do^ i0 , the one loop real- virtual contribution, do^ LO and the double-real 
tree level contribution, dd^ LO - Each component contributes to the same order of a s and so the 
terms with a lower number of loops have a higher number of final state particles thus requiring 
integration over a larger dimensional phase space. The dijet cross section can be written in terms 
of these quantities and expressed as [6], 
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where the requirement that exactly two jets are observed is implemented by the jet function, J% 
which is defined to construct two jets from n final state partons, 
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(1.3) 



The calculation of each of the three terms in Eq. ( \l.2\ > presents different challenges: The 
double-virtual contribution has only explicit poles arising directly from the loop integrations but 
the loop integration itself is highly non-trivial. On the other hand the matrix elements for the 
double-real contribution have been known for many years but their intricate implicit IR singularity 
structure and more complicated phase space integration has caused a bottleneck in the calculation. 
The real-virtual contribution contains both implicit and explicit IR singularities. In this talk we 
will focus on the double-real radiative corrections to the dijet cross section for processes involving 
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quarks and gluons in either the initial or final state. IR divergences are inherent to the double-real 
cross section and arise from field configurations where partons become soft or collinear in any 
combination. The factorization of the colour ordered matrix elements in these IR configurations is 
universal [7] and this is exploited to isolate the singular behaviour of the cross section. 

2. Antenna Subtraction 



A significant effort has been applied to the calculation of the NNLO dijet cross section [8-14] 
and all necessary matrix elements are now known. The remaining challenge is to isolate the im- 
plicit (local) IR singularities in the double-real and real-virtual contributions, to explicitly cancel 
the IR poles between terms and to numerically integrate the resultant finite integrand over the ap- 
propriate phase space. Many techniques have been developed for the purpose of isolating IR poles 
both at NLO and NNLO [15-19]. The antenna subtraction method [18] is one such formalism and 
has been applied successfully at NNLO to the calculation of e + e~ — > 3j [20,21]. The method con- 
sists of constructing a subtraction term, da^ NLO , from antenna functions and reduced multiplicity 
matrix elements which mimics the IR behaviour of the physical matrix elements in all singular 
configurations. A subtraction term is also constructed to absorb the implicit IR singularities of the 
one loop real-virtual contribution. With the subtraction terms defined properly and including mass 
factorization terms, Ao^Iq, the cross section may be rewritten in the form [6], 
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In doing so, the difference between the double-real emission matrix elements and the subtraction 
term is rendered IR finite and the integrand can be integrated over the four-parton phase space, 
yielding a finite result. The subtraction term is constructed in such a way to allow for analytic 
integration of the antenna function containing all IR divergences over the antenna phase space, 
a subspace depending only on the partem momenta involved in a singular configuration. Upon 
integration in d = 4 — 2e dimensions the IR poles may be extracted, cancelled against the explicit 
poles coming from the virtual amplitudes and the finite remainder numerically integrated. The final 
goal is to rewrite the cross section in the form, 
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where each term is individually finite and can be numerically integrated over the appropriate phase 
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space in four dimensions, 

aa NNLO — a<J NNLO ~ NNLO ~ / NNLO 

■Ul 

d&NNLO = ~~ do NN l — / do^ i0 — / da NNL0 (2.3) 

where Xn denotes the phase space corresponding to the parton momenta involved in a singular 
configuration with n unresolved partons. 

In the antenna formalism, the subtraction term is constructed from an antenna function con- 
taining all the IR singular behaviour and a finite reduced matrix element. The antenna function's 
arguments are the momenta of two hard radiators and one or two momenta that are allowed to be- 
come unresolved. The arguments of the matrix element used for the subtraction term are obtained 
from the original momenta via a phase space map [11,22]. The momenta involved in the antenna 
function are mapped onto two composite momenta such that the matrix element is independent of 
the unresolved momenta in the singular limits. The subtraction term for a double real emission n 
parton process can be schematically represented as a sum over terms of the form, 

d&NNLO ~ Xm+lik , ■ ■ ■ , im+l) <8> ^n-mUl > - ' - ) jn-m) (2-4) 

where X m+ 2 is a (m+2) parton antenna function containing at most m unresolved partons and M%_ m 
is a squared matrix element summed over colours and helicities. At NNLO we are interested in 
three and four parton antenna functions to absorb the single and double unresolved divergences. We 
consider colour ordered amplitudes because of their universal factorization properties in unresolved 
limits; the subtraction term is constructed to capture this factorized behaviour. 

F I 

(b) Y 6 (c) Z 6 

Fi gure 1 ! The three hexagon topologies for 2 — )• 4 scattering. / and F denote initial and final state edges respectively. 

A novel way to understand how to construct subtraction terms is to consider the problem dia- 
grammatically. Consider a rc-parton squared matrix element as a n-gon with each edge representing 
the momenta of an external leg. For six parton scattering the real matrix elements are represented 
as hexagons as in Fig. [[]. Two of the six edges are labelled with an / to denote that these edges 
represent initial state partons. There are three independent ways two /'s can be distributed about 
the hexagon, defining three kinematic topologies, X(,,Y(,,Z(„ see Fig. [jj 

As well as the six parton matrix elements, the calculation involves four and five parton matrix 
elements in the subtraction terms. These are represented as box and pentagon diagrams and dis- 
played in Fig. |2[ In these cases there are only two independent kinematic topologies, Xl,F4 and 
X$ , Y$ for the boxes and pentagons respectively. 

The antenna functions which reflect the IR behaviour of squared matrix elements and are 
colour ordered can also be represented in pictorially. Three and four parton antennae are repre- 
sented as triangles and boxes where now the number of initial state edges can be zero, one or two, 
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Figure 2: The two box and pentagon topologies representing four and five parton matrix elements. 

depending on whether we are representing a final-final, initial-final or initial-initial antenna, see 
Fig- 1 
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Figure 3: The three triangle and two additional box topologies used to represent the three and four parton antenna 
functions. 

Using this diagrammatic approach, the factorization of the mauix elements is quite intuitive 
and follows a diagrammatic algorithm, represented in Fig. [| First consider a final state parton be- 
coming unresolved and identify its hard colour adjacent neighbours. The hard neighbouring edges 
are then stretched and pinched to a point, separating the original polygon into two daughter poly- 
gons. The daughter polygon containing the unresolved edges represents the appropriate antenna 
function and the remaining polygon containing only resolved edges represents the reduced matrix 
element with mapped momenta. 





a • 'ty 



Figure 4: A diagrammatic representation of single unresolved parton factorization for six parton scattering belonging 
to the X(, topology. The unresolved edge (blue) is adjacent to two hard edges (red). Edges are labelled / and F to denote 
initial or final state partons. 



The unresolved limits of the six parton matrix elements can be systematically examined by 
simply taking all the possible pinches of the appropriate hexagon. In the single unresolved limits 
the hexagon is pinched to give a product of a triangle and a pentagon, see Fig. A colour 
connected double unresolved limit is obtained by letting two adjacent edges go unresolved; in this 
case the hexagon is pinched to form the product of two boxes. Almost colour connected limits are 
derived by iteratively pinching out two edges separated by a single hard edge, yielding a product 
of two triangles with one edge in common and a box. Colour disconnected double unresolved 
limits are taken by pinching out two edges with no neighbours in common resulting in the product 
two triangles with distinct edges and a box. The four parton antenna functions also contain single 
unresolved limits of their own which have to be removed. This is achieved by applying the same 
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algorithm to the box as was done for the hexagon; single unresolved edges are pinched out and 
subtraction terms are constructed from a product of two triangles. In this way the full subtraction 
term can be constructed from blocks of subtraction terms which remove unresolved limits of the 
matrix elements in a systematic way whilst not introducing spurious singularities of their own. In 
Fig. |5| a block is defined to remove a specific double unresolved limit from the six parton matrix 
element using a four parton antenna function whose spurious single unresolved limits are removed 
in turn by products of triangles. 
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Figure 5 : The block of subtraction terms used to absorb divergences arising from colour connected partons F\ and Fi 
(blue) becoming unresolved in a six parton matrix element (hexagon) belonging to the X(, topology. Colour connected 
double unresolved limits are removed using the four parton antenna (box). Spurious single unresolved limits of the 
box are removed using iterated three parton antennae (triangles) including the extra limit when F3 becomes unresolved 
(green). 



The specific antenna functions and matrix elements used in this block depend on the par- 
tonic channel under examination and can be read off directly from the diagrams when the edges 
are labelled with the appropriate parton species. In six gluon scattering, the block of subtraction 
terms in Fig. |5] absorbs divergences coming from the i,j unresolved limit of the matrix element 
M®(l,2,i,j,k,l), where the hat denotes initial state partons 1 . Translating the block of diagrams 
from Fig. || into the specific functions for this process yields the subtraction term [6], 

F°(2,i,j,k)M° 4 (\Xm,l) - J§&iJ) F^2M,k)M^\X(m),l) 

- j§(i,j,k) F^2,Qf), {jk))Ml(\X{{mfk))J) 

- f 3 {j,k,2) F°(%i,(jk))M°(iX(iUQ),l) (2-5) 

where f^,F®,F® are the various three and four gluon antenna functions, and their arguments reflect 
the parton momenta generated via the appropriate phase space maps [11,22]. This process was 
considered in [6] and such a block of subtraction terms can be found in the appropriate double-real 
subtraction term formula. The generality of this approach can be seen by considering a different 
physical process, for example a quark-antiquark pair scattering into four gluons where we consider 
the colour ordered six parton matrix elements M®(l q ,i,j,k,l,2q). In this case the block of diagrams 

1 This set of matrix elements is represented diagrammatically by a hexagon belonging to the X b topology 
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in Fig. |5| translates to a subtraction term given by, 

Dl{l q ,i,j,k) M2(l q ,(ijk),l,2 q ) - 4(\ q ,iJ)D%\ q M^)Mt(\ q ,{^),l 1 2 ij ) 

-f$(i,j,k) D° 3 (l q M,W)M^\ q ,((i^k))J,2 q ) 

- ^(l 4> *,jO^(i«,i,(^)^(I«,(W),'i2 ? ) (2.6) 

where d®,D®,D® are quark-gluon antenna functions defined in [18]. By deriving blocks of subtrac- 
tion terms in this manner and carefully keeping track of any over-subtraction, the full subtraction 
term can be formulated systematically. The details of the subtraction term for each partonic chan- 
nel vary according to the matrix elements and antenna functions involved but can be derived within 
this basic framework. In order to keep track of the parton species involved in a process it is often 
convenient to draw the polygon edges to reflect the parton it represents as in Fig. |(| in this way the 
antenna species and specific singular limits can be read off the diagrams directly. 




k k 



(a) (b) (c) 

Figure 6: Diagrams can also be drawn with wavy edges to represent gluons and straight edges to represent quarks for 
either matrix elements or antenna functions, (a) M^{\ q ,i, j ,k,l,2q), (b) F®(i,j,k), (c) Dfy(i q ,j,k,l) 

3. Conclusions 

In this talk, I have outlined the general procedure by which single and double unresolved IR 
singularities can be isolated using the antenna subtraction method. The method relies on two fac- 
torization theorems: antenna phase space factorization and universal matrix element factorization. 
These properties are exploited to construct a subtraction term which mimics the double-real ra- 
diation cross section in the singular limits and thus absorbs all divergences without introducing 
additional spurious singular behaviour. A novel pictorial method is introduced for constructing 
the double-real radiation subtraction term within the antenna subtraction formalism. The method 
allows general blocks of subtraction terms to be generated using only kinematic information; these 
can then be dressed with process specific information to determine the final form of the subtraction 
formula. 
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